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Abstract. Darboux Wronskian formulas allow to construct Darboux trans- 
formations, but Laplace transformations, which are Darboux transformations 
of order one cannot be represented this way. It has been a long standing 
problem on what are other exceptions. In our previous work we proved that 
among transformations of total order one there are no other exceptions. Here 
we prove that for transformations of total order two there are no exceptions 
at all. We also obtain a simple explicit invariant description of all possible 
Darboux Transformations of total order two. 
Mathematics Subject Classification Numbers: 53Z05, 35Q99. 



1. Introduction 

Classical Darboux transformations and their generalizations are methods for 
obtaining analytic solutions of linear Partial Differential Equations (PDEs) . They 
also serv e as a leverage for l arger theories for solution of non-linear PDEs, see for 



example iMatveev and Salld (|l99lh and references therein. 

In the present paper we are concerned with the intertwining relations J\f o C = 
L\ o A4 for operators of the form 

(1.1) L = D^Dy + aD^ + bDy + c , 

where the coefficient s may be non c onstant. Since they have been introduced in 
the classical work of Darboux ( 1889I ). we shall call them Darboux transformations 



too. PDEs corresponding to such opera t ors ap pear also as part of the problem of 



the search of flat metrics, see iKricheveii (|199 



r 



Given Linear Partial Differential Operator (LPDO) C and some LPDO M., the 
coefficients of the resulting operator Ci and of the auxiliary operator M can be found 
algebraically. There are two choices of A4, which always lead to a DT for a given 
operator (jl.ip : — + b, and Ai — Dy + a. These Darboux transformations 
have a special name; Laplace transformations. The latter are not be confused with 
Laplace integral transforms. 

There is also a large class of Darboux transformations generated by operators M 
that are constructed using so-called Darboux Wronskian formulas. These are based 
on the assumption that we know some number of linearly independent particular 
solutions of the initial PDE, Cij] = 0. This class is a very large class and Darboux 
transformations of arbitrary orders can be constructed provided we know enough 
number of particular solutions. Laplace transformations, which are Darboux trans- 
formations of order one do not belong to this class. 
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Lap lace transformations are particularly good Darboux transformations, see lTsarev 



(|2005[ ). and they have been the only known examples of Darb oux transformations 



that cannot be described by Darboux Wronskian formulas. In IShemvakoval (j2012r ) 
we have proved that a Darboux transformation of total order one is either described 
by Darboux Wronskians or is a Laplace transformation. The problem reduces to so- 
lution of a non-linear PDE. The PDE was not so large and noticing some interesting 
structure we were able to tackle the problem. 

After that work it was still unclear whether there are some exceptional trans- 
formations, that is such that cannot be described by Darboux Wronskian formulas 
among Darboux transformations of orders higher than one. This problem is reduc- 
ing to solution of a system o f two large non-linear PDEs, for which methods of the 



previous work IShemyakoval (j2012[ ) were hard to apply. We, however, succeeded in 
proving that Darboux Wronskian formulas complete for transformations of order 
two in a different and rather elegant fashion, and present this proof in this paper. 

Recently there have been s everal n ew ideas to tackle Darboux transformations 
and related problems. Thus, iTsarev and Cassidv and S inger (2010) have 



made very important progress in the description of factorizable operators corre- 
sponding to linear PDEs in terms of certain abelian categories and algebraic groups, 
respectively. In the present paper, we adopt an approach that is based on ideas of 
Differential Geometry, and is constructive. 

Our main result is an elegant proof that all Darboux transformations of total 
order two can be described by Wronskian formulae (Theorem 17. 5p . The second 
achievement is an easy to use invariant description of all these Darboux transfor- 
mations (Theorem 17. 4|) . 

The paper is organized as follows. Darboux transformations of total order two 
are those that has in one of the following forms: 

M = m2oDxx + mioDj; + moo , 
M = mo2Dyy + moiDy + mQo , 
M = mioD^ + moiDy + mQQ , 

where the rriij € K are not necessarily constant. In Sec. [3] we show that to cover 
all Darboux transformations of total order two it is enough to consider M. is the 
form Dx + qDy + r, where p and q are some functions. After some preparation 
in the next two sections, we introduce in Sec. |6] new transformations of the pair 
{C,M}, which we name gauged evolution. We determine the generating invariants 
uniquely defining the equivalence classes under these transformations and use them 
to invariantize the nonlinear system of PDEs defining all possible Darboux trans- 
formations of total order two. The invariantized system is easier and can be solved 
explicitly by classical methods, however, even though we have a technical solution, 
it is in quadratures and it is useful neither for invariant description of Darboux 
transformations, nor to judge whether Wronskian Formulae give all such Darboux 
transformations or not. 

Therefore, we need a further invention, Theorem 17.31 through which we are able 
to obtain an elegant general solution (Theorem 17.41) of the invariantized system of 
PDE. We still have to remember that even if the invariantized system of PDEs has 
solutions, the existence of Darboux transformations depends also on the existence of 
a solution of a nonlinear PDE system (|6.3p . where we return from gauged evolution 
invariants to the coefficients of operators C and Ai. In the proof of Theorem [73] we 
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resolve this problem and conclude that for every Darboux transformation of total 
order 2 there exist two linearly independent partial solutions of Cu = 0, such that 
it can be constructed using Darboux Wronskian formulas. 



2. Preliminaries 

Let K he a. differential field of characteristic zero with commuting derivations 
dx: dy. Let K[D] = K[Dx,Dy] be the corresponding ring of linear partial differential 
operators over if, where Dx,Dy correspond to derivations dx,dy. 

Operators C G K[D] have the general form C = J2t+j=o ^^ij^x^i: where aij G 
K. The formal polynomial Sym£ = aijX^Y^ in some formal variables X, Y 

is called the symbol of C. 

One can either assume field K to be either differentially closed, in other words 
containing all the solutions of, in general nonlinear. Partial Differential Equations 
(PDEs) with coefficients in K, or simply assume that K contains the solutions of 
those PDEs that we encounter on the way. 

Definition 2.1. An operator Ci G if[-D] is called a Darboux transformation of an 
operator C G ^^[f], if Sym(£) = Sym(£i), and there exist operators JV E K[D] 
and M e K[D] such that 

(2.1) AfoC^CioM. 

In this case we say that this Darboux transformation corresponds to pair {C,A4}, 
and that operator Ci is associated with, or Darboux- conjugated to operator L and 
use the notation 

Ci = ip{C,M,N) ■ 
Note that coefficients of the operators are not required to be constants. 

Darboux transformation implies the following transformations of kernels: Ker C ^• 
Ker>Ci : V ^ M{ip) and requires Sym(7V{) — Sym(A/'). 

Definition 2.2. The Darboux transformation of an operator (|l.ip . where a,b,c are 
not required to be constants, is called a Laplace transformation if the corresponding 
operator A4 is either A4 = Dx + 6, or = Dy + a. 

Laplace transformations are the most well-st udied cas e of Darboux transforma- 
tion and have several important prop erties, see iDarbouxl ([1889). 

One of the most famous results in iDarbouxl |T889[) concerns Darboux transfor- 
mations for operators of the form (jl.ip and can be formulated as follows. 



Theorem 2.3 (Darboux). Let we C be an operator of the form (jl.ip andtpi, 
Kcr£ be linearly independent then 

(2.2) A^(V') = W^™,„(V',V'l,...,V'm+n) 

defines some Darboux transformation for operator C. 

Here Wm,n is a Wronskian-likc function Wm,n{'4', i^i, ■ ■ ■ , ipm+n) — 



m-\-7i 



Dx^ 
Dxil^i 



Dyi, 
Dylpi 



Dyip 



y Vm+n 
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That is a Darboux transformation of order m + n ca n be built using m + n par- 
ticular solutions of the initial equation C{ip) — 0, see IParbouxl ( 18891 ). Darboux 



Wronskian-like formulas (|2.2p provides a large class of possible Darboux transfor- 
mations. 

3. Normalization of Darboux transformations 

3.1. Normalization of Darboux transformations using expansion. 

Lemma 3.1. Let C be of the form (jl.ip and M. he an operator of arbitrary order 
from K[D\. Let M. define at least one Darboux transformation for operator C, then 
for any given operator A G ^^i^*] there exists also a Darboux transformation for 
the same C with M = Ai + Ao C. 



Proof. Equality (|2.1|) implies that 

£io{M + AoC)^Mio£ + CioAoC^ {Mi + CioA)oC 

is true for an arbitrary operator A G -^''[i^]- Therefore, there exists a Darboux 
transformation for C with Ai = A4 + A o C. □ 

Definition 3.2. Lemma lOl describes transformations of pairs of operators {C, M}. 
It shows that such transformations preserve the property of the existence of Dar- 
boux transformations for a given operator C, and splits the operators Ai into 
equivalence classes. We name this transformation an expansion. 

Remark 3.3. Notice that the resulting operators of the initial Darboux transfor- 
mation and of the one generated ioT C hy A4 + Ao C are the same. 

Given operator C of the form ()1.1|) . we shall be considering different pairs 
{C,A4}, where A4 are operators in K[D]. Using expansion we can eliminate all the 
mixed derivatives in M in the case of C of the form (|l.ip . 

Definition 3.4. Let C G K[D] be of the form (|l.ip and M G K[D] be an arbitrary 
operator, then denote the result of elimination of the mixed derivatives in Ai using 
£ as TTciAt). 

Definition 3.5. Given C G K[D] of the form (II. ip . we define the bi-degree 
deg£ Ai = (to, n) of operator Ai with respect to C as follows: m is the highest 
derivative with respect to in nciAd) and n is that with respect to By. We shall 
say that m + n is the total degree of Ai. 

Definition 3.6. By the degree or total degree of a Darboux transformation Ci ~ 
(p{C,Ai,N') we shall understand the degree or the total degree of Ai. 

3.2. Normalization of Darboux transformations using composition with 
Laplace transformations. 

Definition 3.7. Let there be a Darboux transformation of arbitrary C G K[D] 
defined by some Ai G K[D], that is (|2.ip holds. Let the result, operator Ci G K[D] 
be transformed into some £2 G K[D] by a Darboux transformation defined by some 
A^i G ^[-D], that is A/i o£i = £2 o A^i for some A/i e K[D]. Then the composition 
of these two Darboux transformations is a Darboux transformation transforming £ 
into £2 defined by A/i o A/" o £ = £2 o A^i o Ai. 

The following lemma allows us to use expansion and composition together. 
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Lemma 3.8 (Correctness of the composition of two Darboux transformation with 
expansion). The result of composition of two Darboux transformations does not 
depend on the choice of the operator A4 within its class of equivalence under ex- 
pansion. 

Proof. Let there be a Darboux transformation of arbitrary C G defined 
by some M. G K[D], i.e., (j2.ip holds. Let the result, operator Ci G K[D\ be 
transformed into some £2 G ^[^] by a Darboux transformation defined by some 
Ml G K[D], i.e., TVi o £1 = C2 oMi for some TVi G K[D]. Consider {M+AoC) 
and M.I + B o Ci for some A,B ^ K[D], which belongs to the same classes of 
equivalence under the expansion as A4 and M.i correspondingly. That is we have: 

{N + CioA)oC ^ Cio{M+AoC), 
{Mi+C2oB)oCi = C2o{Mi+BoCi) . 

Then the composition is 

{Ml + £2 ° B) o {JV + £1 o A) o £ ^ {Ml + C2 o B) o £1 o {M + Ao C) , 

which using equality p.ip can be re-written as 

{Ml + £2 o B) o {M + £1 o A) o £ ^ £2 o {Ml + B o £1) o {M + Ao £) . 

After expanding some multiples and re-grouping we have 

{MioM + £2°C+Mio£ioA)o£ = £20 {MioM+£o£ + Bo£i o M) , 

whereC = BoM+Bo£ioA, and f = MioA+Bo£ioA. Finally, substituting A/'o£ 
instead oi £ioM we obtain that the "TW" operator of this Darboux transformation 
belongs to the same equivalence class that Mi o M does under the expansion 
transformation. □ 



Then one of the results of lDarbo"ux ( 1889h can be interpreted as follows 



Theorem 3.9. Let £ G K[D] he of the form and M G K[D] define a Darboux 

transformation for £, and deg^ M = {m, n). Let Mx — Dy + a and My = + b 
define LTs for operator iLl]) . Then 

(1) AegT:{M o Mx) = (m - 1), 

(2) degi:{MoMx) = (m + l,n- 1), 

(3) tt{Mx ° My) — by ~ c + ab, which is an operator of order zero, 

(4) T:{My o Mx) = o,x — c -\- ab, which is an operator of order zero. 

Summarizing all the results we can formulate the following theorem. 

Theorem 3.10. Let £ G K[D] be of the form I7j]) and M G K[D] define 
a Darboux transformation for £. Let deg£ ~ {m,n). Then for every i = 
1, . . . ,min(rn,n), there exists an operator Mi without mixed derivatives having the 
property that deg£ Mi = {m — i,n + i) and Mi defines a Darboux transformation 
of£. 

Lemma 3.11 {M can be multiplied by a function on the left). Let there exist 
a Darboux transformation of operator £ G i^ii*] with some operator M G K[D]. 
Then for every invertible element p K there exists a Darboux transformation of 
operator £ with operator pM . 
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Proof. The conditions of the lemma imply that, for some J\f,Ci G K[D], equal- 
ity (|2.ip holds. Therefore, p o J\f o C = p o Ci o p^^ o p o is true also. Since the 
symbol of £i is not altered under gauge transformations, then operator po/3x °p~^ 
is an operator of the form (jl.ip and we have proved the statement of the lemma. □ 

Let £ € K[D] be of the form and A4 £ K[D] of arbitrary form and order 
defining a Darboux transformation for C Theorem l3.10l and Lemma fS .111 imply that 
using operations of expansion, composition with LTs, and division by a function on 
the left, we can bring such Darboux transformation into a normalized form with 
M having no mixed derivatives and having one of the following symbols: 

Sym(A^) ^ fc > , 

Syni(A^) = y^ fc > , 
Sym(yW)= X'= + gy^ fc > 0, g 7^ . 

Before we decide which of these to use in further considerations, let us consider the 
uniqueness problem for Darboux transformations. 

3.3. Uniqueness of Darboux transformations for given C and Ai. 

Theorem 3.12. Let C e K[D] be of the form U.l\) and let M G K[D] define some 
Darboux transformation. Then, unless for its normalized form we have Sym(A^) = 
or Sym(A^) = , such a Darboux transformation is unique. 
If for its normalized form Sym(A^) ~ (corresp. Sym{M) — ) and there 
exist two Darboux transformations: Ci — ip{C, Ai^Af) and Ci + £[ = ip{C,Ai + 
M',^f + ^^'), then 

Sym{M[)^X'^-\L[^Dy+j , \ 
^ ' {corresp. Sym{M[) = Y''-^,L[ = D^+-f) . j 

for some j G K . 

Proof. Since C is of the form (|l.ip . for C[ there are only four possibilities: 

£[ = D, + l3Dy+j, /3 7^ , 

£[= 1. 

Ci = (p{C,M,JV) and Ci + C[ = (p{C, M + M' ,M + M') implies 
(3.2) M[oC^C[oM . 

Case C[ = 1 cannot take place because if it does then j\4[ o C = Ai, which is 
impossible as Sym(M) cannot be divisible by XY. 
Let Sym(X) ^X^,k> 0, then 

Sym(7U;) -XY ^ Sym(/:;) • X^ , 

which implies that Sym(£']^) must be divisible by Y , which is only possible if C'l = 
Dy + 7. Then Sym(A^']^) cannot contain extra Y-s, and therefore, Sym(A^'j^) = 

x^-\ 

Analogously, if Syni(A^) = F'"', fc > 0, we have C[ ^ + ^ and Sym(A4i) = 

yfe-l 

Let Sym(M) ^ X'' + qY^ , then 

Sym(X'i) -X-Y^ Sym{C[) ■ {X^ + qY'') , 
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which means that S'ym{C'i) must be divisible by XY, which is impossible. □ 

Theorem 13.121 guarantees uniqueness of a Darboux transformation for given A4 
and C if Sym(X) = X'' + qY^ . Further below we shall be interested in Darboux 
transformation of the total degree two, and we choose the normal form for such 
transformations with 

Sym(X) ^X + qY . 

4. Existence of a Darboux transformation defined by M of 

bi-degree (1, 1) 

Even for the simplest case of M. of total degree 1, the problem of describing 



all Darboux transformations is not easy IShemvakoval (j2012l) . For the case of M. of 



total degree 2, which is considered here the problem becomes very difficult. 
Theorem 4.1. Let C he of the form il.l]] and M G K[D] in the form 
(4.1) M = Dx+qDy+r . 

If there exists a corresponding Darboux transformation, then the corresponding op- 
erator J\f is given by Af — Ai ~ {In q)x + qy ■ The necessary and sufficient conditions 
for the existence of a Darboux transformation for such pair (£, A^) are 



-qrx + q^ry + qxr - bqx + bxq + q'^{by - aqy - Ox) ~ q^ay+ 

+qyqx - qxyq = 

-cqx + (c - ar)qyq + {ar + ry)qx + {cy - ray)q^+ 

+ {rry - avx - ryb - rxy - rax + Cx)q = 



Proof. Compare the corresponding coefficients on the both sides of equality (I2.1|l . 

□ 



Darboux theorem 12.31 provides us with a particular solution of the system (|4.2p . 
The following statement is Theorem 12.31 written out more explicitly for the case of 
M of bi-degree (1,1). 

Theorem 4.2 (Darboux main theorem for bi-degree (1,1)). Let C E K[D] be an 
arbitrary operator of the form (jl.ip and V'l, V'2 be two linearly independent solutions 
of Cip = 0. Then there exists a Darboux transformation with 

a B 
M^Dx + - Dy + ^ , 
a d 

where 

d = -i>l'4>2y + ^2'4'ly , 

a = lpl'lp2x - i>2i>ix , 

P = ~'4'2x'4'ly + ll^2y1plx ■ 

Remark 4.3. If we denote by ■0 the ratio of these particular solutions, 

V' = — , 

then A4 in the statement of Theorem 14.21 can be written in more simple form: 

AA - n _ n J- ^ij^^^ _ 

■' ' — I ^y ' ) I I ■ 

Wy WlVy Vl 
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Remark 4.4. In order to describe all Darboux transformation for of bi-degree 
(1, 1) we need to solve system (|4.2p for g,r, where a,b,c are known and are not 
constants in general. Usual differential elimination techniques does not lead to a 
general solution. 

A different approach can be to notice that in the system (|4.2p the second equation 
is non-linear in both q, r, while the first equation is nonlinear in q only. The first 
equation is a linear first-order non-homogeneous PDE on r, and since we know its 
particular solutions (Theorem 14. 21) . one may solve it in quadratures. These quadra- 
tures are expressed in terms of q, and therefore, after substituting the expression 
for r into the second equation one gets even more nonlinear, rather large, PDE. 

In the rest of the paper we shall be proving that the general solution of sys- 
tem (|4.2p is given by the class of particular solutions from Theorem 14.21 

5. Gauge Transformations of Pairs and Corresponding Invariants 

Our plan is to address our problem using invariants methods. In this section 
we study gauge transformations of pairs which are almost classical with 

the only difference that we apply them to the pairs of operators. These transfor- 
mations are not strong enough to simplify our system significantly, and completely 
new transformations will be introduced in Sec. |6l However, we shall use gauge 
transformations of pairs too. 

Definition 5.1. Given some operator TZ £ K[D] and invertible function g ^ K, 
the corresponding gauge transformation is defined as 

where o denotes the operation of the composition of operators in -?ir[r']. It is 
convenient to take g in the form g = exp(a). Then we shall avoid fractions while 
writing this transformation out on the coefficients of TZ. 

Our first step towards simplification of the problem is the following simple ob- 
servation. 

Lemma 5.2. Let C = D^Dy+aD^+bDy+c € K[D] and M = D^+qDy+r € K[D]. 
If a Darboux transformation exists for the pair (C,A4), then one exists also for 
{A4^,C^), where g is an arbitrary invertible element of K . 

Proof. Indeed, from the Darboux equality (|2.1|) for the pair (A^, £), we have 

g^^ ojVogo o Co g — g^^ o Ci o g o g^^ o M. o g , 

and, therefore, o £f = £f o M.^ . Recalling that gauge transformations do not 
change the symbol of an operator, we conclude the proof of the lemma. □ 

Therefore, it is natural to consider our problem for the equivalence classes of the 
pairs (A4, £). In order to define every class uniquely we determine a generating set 
of all the invariants of these pairs under the gauge transformations. 

Definition 5.3. Let TZ G K[D] be an operator and T be some transformation 
acting on Then a function of the coefficients of TZ and of the derivatives of 

these coefficients is called a differential invariant if it is unaltered under the action 
of T on TZ. 

The sum, and the product of two differential invariants is an invariant, as well as 
a derivative of an invariant is also an invariant. In the infinite set of all possible 
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differential invariants there is some subset (not necessarily proper) of differential 
invariants which generate all others using algebraic operations and derivatives. Such 
a subset we shall call a generating set of invariants. 

Theorem 5.4. Let C D^Dy + aD^ + bDy + c e K[D] and M = + qDy + 
r e ^[-D]. On the set of all pairs (£, A^) of such operators consider the gauge 
transformation of those with function exp(a) ; 

The following functions are invariants and in addition form a generating set of all 
differential invariants for such transformations: 

(5.1) m = a-by, \ 

n = ab — c + Ox , 
R ~ r — b — qa . 

Remark 5.5. Functions h = ab — c + ax, k = ab — c + by arc known as h- and 
/c-Laplace invariants as they are invariants of operator C considered individually 
(without Ai) under the gauge transformations. Both of them are present here: h 
is present in its original form, and k is hidden in m as m — h — k. 

Proof. To find a generating set of differential invariants we use the method of reg- 
ularized moving frames introduced by Fels and Olver in Pels and Olveil (Il998l). A 
good overview of recent developments in the area can be found in lMansfieldT i 2010[ ) . 
Note that our case is infinite dimension al, so the connected difficulties have been 
treated in I Olver and Pohianpelto ( 20051 ). 

The transformations in question can be defined coordinate-wise as follows. 

ai = a + Uy, 
bi ^b + ax, 

ci = c + aux + bay + Uxy + axUy , 
n = r + ax + qay , 

where £<='<^p(") = DxDy + aiDx + biDy -f ci and A4<=''p(") = Dx + qiDy + n. We 
are choosing a cross-section as follows 

(5 2) 

where J is a string of the form x . ^.x^y . . . y, where n — 0, 1, 2, ... , and where 

m = 0, 1, 2, . . . and X is a string of the form x . ^. x, where I — 1,2, ... . For every 

I 

f G J'^li'], the notation fj stands for the mixed derivative of /: of order n order 
with respect to x and m with respect to y, while fx stands for the lih derivative 
of / with respect to x. 
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This gives us non-contradictory all the values for the parameters of the pseudo- 
group action, ax,ay,axx, ■ ■ ■■ 

ax = -b , 

^xy — O^x 5 

while the value for a we choose arbitrary, as it does not appear explicitly in the 
definition of the pseudo-group action. Then we evaluate the edge invariants on the 
frame: 

{bi)y = by + Uxy = by - ax , 

ci ~ c ^ ab — ab — Ux + ab = c — Gx — ab , 

ri = r — b — qa , 

which constitute the generating set of differential invariants of the pair under the 
gauge-transformations of the pair. Invariants m and h differ by a sign from the 
first two we have just obtained, the third invariant we have obtained is exactly R 
from the statement of the theorem. □ 

Definition 5.6. We shall call invariants (15. ip the gauge invariants of the pair. 

Express the coefficients of the pair {C,M) in terms of invariants: 

r ^ b -\- qa + R , c = ab — h + ax ■ 



After these substitutions system (|4.2|) does not depend on b itself, but only on 
its derivatives by and bxy Therefore, we can effectively use gauge-invariant m by 
enforcing substitution 

by = ax — m . 
Then system (|4.2|) simplifies to the following one. 

fla + qxh - qhx - q^hy - qxm + qxRy+ 

+ qrux — qRxy — qyqh — qRm + qRRy — 

where ft — — 2g^m + q^Ry + qxR + qyqx — qRx — qxyq- Therefore, this system can 
be simplified further: 

- 2q'^m + q^Ry + qxR + qyqx - qRx - qxyq = 



(5.3) qxh — qhx ~ q hy ~ qx^n + qxR 



y 



+ qnix — qRxy ~ qyqh — qRm + qRRy = , 

In our problem Darboux transformation, operator C is considered to be given, 
therefore, the gauge-invariants h and m are given, and the problem is reduced to 
the search of the general solution for system (15. 3p with respect to unknowns q and 
R. 

Although system (|5.3I) is visually shorter than system (|4.2I) . it is still hard to 
solve using the usual methods, such as the differentiation-cancellation technique. 
Note that invariantiation and in particular moving frames method hav e been 



useful for investigation of Darboux like me t hods earlier, see for example lOlvei 
(|2nilh . lVassiliOT] (|2nilh . lGomez-Ullate et al\ (|2nilh . 
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6. Gauged Evolution of Pairs and Corresponding Invariants 

Lemma 6.1. Let C,A4 E K[D] be two arbitrary operators in K[D] and let C be 
of order larger than one, while M. is a first-order operator. Then if a Darboux 
transformation exists for the pair (£, Ai), it also exists for the pair (£ + PAd, Ai), 
where /3 € K is arbitrary. 

Proof. The existence of a Darboux transformation for the pair (£, Ai) means that 
for some N,Ci G if [-D], where £i has the same symbol as £. Therefore, 

N oC^ CioAd . 



Then 



(6.1) 



Mo{C + (3oM)^CioM+Mol3oM = {Ci+Mol3)oAi. 



Since M must be of the same order as A^, TV is a first-order operator. In addition, /3 
is zero-order operator. Therefore, the symbol of the operator Ci +Afo f3 is the same 
as the symbol of £i, which is the same as the symbol of C. Therefore, equality (|6.ip 
defines a Darboux transformation for pair {C + (3Ai, Ai). □ 



Definition 6.2. These transformation on the pairs, that is {£,Ai) — {£,Ai): 

C C + l3oM , 

M -f M . 

we shall call evolution of the pair (or ^-evolution of the pair). 

Definition 6.3. Let C = D^Dy+aD^+bDy+c € K[D] and Ai G K[D] is arbitrary. 
On the set of all the pairs {C,Ai) of such operators consider the consequential 
application of the gauge transformations and of the evolution: for given a, (3 E K: 



(6.2) 



We shall call these transformations gauged evolution of the pairs. 

Theorem 6.4. Let C = D^Dy + aD^ + bDy + c e K[D] and Ai = + qDy + 

r e K[D]. The gauged evolutions of the pairs {C,A4) of such operators have the 
following generating set of differential invariants: 



(6.3) 



h = q , 
h = 2m - R 



h = 2h 



IV 

1. 



Remark 6.5. Notice that gauged evolution generating invariants (|6.3p are ex- 
pressed in terms of generating gauge invariants q, h, m, R only. This means that 
the gauged evolutions split the set of pairs {C,Ai) into larger equivalence classes 
than the gauge transformations of pairs do. Also we see that those "small" gauge 
classes can belong to the "larger" gauged evolution classes only entirely. 
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Proof. Evolution (|6.2p can be defined coordinate-wise as follows: 

ai = a + ay + P , 
h ^ b + ax + , 

Ci = c + OQfj. + bay + a^y + aya^ + f3r + fSa^ + fSqay , 
91 = 9 , 

where £°'^P(")+/3>1<=''P(") = D^Dy+aiD^+biDy+ci and TW^^p^") = D^+giDj^+ri. 
We are setting a cross-section by setting most of the coordinate functions to zero: 

(ai)j = 
(6i),7 = 

where J and X are the same notations as in (|5.2p . Then at the beginning we have 
three equations, 

ai = , 5i = , ri = 
and three variables, parameters to determine: 

(5 , ax , ay . 

The determinant is not 0, so there is a unique solution for such a system. At the 
next step we consider first prolongations only, which gives us 5 equations for 5 
variables and this linear system has non-zero determinant. In general, considering 
i-th prolongation we have 2i + S variables and the same number of equations, and 
a non-zero determinant of the corresponding linear system. Therefore, we have 
defined a frame, and the generating set of invariants in this case consists of the 
corner invariants: 

rr Tx+bx , qxT - qxb , 

ly = ry ~ Gx 1 qya - qay + by , 

Ox br qxb qxT ab + ar qa? b^ r^, bx 

r = b + qa + R , 
c ~ ab — h + ax , 
ax = by + m . 

we obtain (up to a sign and a multiplication by 2) the invariants claimed in the 
statement of the theorem. □ 

Definition 6.6. We shall refer to invariants (|6.3p as gauged evolution invariants. 

7. Solution of the PDE System. Description of All Darboux 
transformations of total order two. 



Substituting 



In Lemma 16.1 1 we showed that the property of the existence of a Darboux trans- 
formation for a pair is invariant under the gauged evolutions. This does not nec- 
essarily mean that there is some explicit invariant form for system (j4.2p . Theorem 
below demonstrates, however, that in this particular case, we can have such explicit 
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invariant form. We also see that the invariantizing system can be written in much 
simpler form than system ()4.2|) . 



Theorem 7.1 (Necessary and sufficient conditions for the existence of a Darboux 
transformation in terms of evolution). Given pair (C,A4), where C = D^Dy + 
aDx + bDy + c S and = D^ + qDy + r G i^[-D], there exists a corresponding 

Darboux transformation if and only if its evolution invariants (q, 12,13) satisfy the 
following two conditions simultaneously: 

(7.1) l2 + Qxv = , 

(7.2) Is^j; + qh^y + {qy - qx / q) H = QxQxy - Qxxy , 
where Q — Inq. 

Proof. Expressing m and h using the second and the third equations of (j6.3|) and 
using by = Ox — m system (|5.3p can be written as in the statement. □ 



We also invariantize the class of particular solutions for system (|4.2p which we 
derived from Darboux Wronskian formulas: 

Theorem 7.2 (Darboux transformations constructed from Wronskians). Let C = 

DxDy + + bDy + c € K[D] and let i/'i, "02 be two linearly independent elements 

of its kernel. Let tp — and A ~ and B = . Then for C there exists 

Wl Wx Vy 

a Darboux transformation such that the evolution invariants of the corresponding 
pair {M., C) are as follows: 

B 

(7.3) h - By - Ax , 

h-~Ax + ^ 

Proof. Compute the values of the gauged evolution invariants (16.31) for A4 con- 
structed using Darboux formulas given in Theorem 14.21 Then we have 

q 

Vy 



-i>x 














'^-^xx'4^xy 


4'xyy 


i^x 






'4^xxy 


'4^xx'4^xy 


Vxy 


i>x 




2'ipxip 



V': 



2 



which can be re-written in a very short form using notations A and B. □ 

Value of I3 from Theorem l7.2l is a particular solution of (|7.2p . a first-order linear 
non- homogeneous PDE on I3. 

Let us solve (|7.2p . One useful idea is to consider q in the form q — — -, where 

Zy 

z is not required to be a ratio of two particular solutions of C{'ip) ~ 0. For this 
q K, invariants /i, I2 can be computed straightforwardly: Ii — q,l2 = —Onq)xy 
Equation ()7.2p is a first order non-homogeneous, and the general solution can 
be obtained as the sum of its particular solution and of the general solution of the 
corresponding homogeneous PDE. As a particular solution we take the expressions 
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from (|7.3|) . Note that for this particular solution z = ip, where ip is some particular 
solution of C4> — 0. 

Using expression for q in terms of z the homogeneous PDE corresponding to (I7.2p 
can be written as I3 ( — ^ ^_^y_ J _ _|_ —J^ ^ = 0. Considering in the form 

= e'-' for suitable J E K[D] and assuming Zy ^ 0, the PDE is equivalent to 
(7.4) T - + ^ = , 



Zx Zy 

the non-homogeneous part of the equation being T = — ^ • After applying 

^y 

the method of characteristics, which we follow more carefully in the proof of the 
next theorem, Theorem 17.41 we choose change of variables — x^rj = z{x,y), and, 
in the new variables, equation (|7.4p has the form 



and, therefore. 



where y — y{^, 77) is the solution of y) — rj for y and F{rii) is an arbitrary function 
of rj. Changing variables back we have 



J = J T{U{tz))z^{ty{^,z))d( + F{z) . 
This means that the homogeneous part of PDE (|7.2I) has the general solution 

/3 - G{z) -expf f T(e, m, (c, 



where G{z) is an arbitrary function of z. Then since we know the particular solution, 
^30 ("0)1 where ip — '02/V'i for two particular solutions of Cu = 0, the general 
description of all Darboux transformations of total order two will be as follows: 



^x 

1 = : 

Zy 

Z.r 



(7-5) '^-[^ 



G(z) -exp (^j T{U{Lz))H{Lm,z))ddj +4o(V'2M) • 



However, formulae (j7.5l) does not serve either of our two purposes: short description 
of all possible Darboux transformations and proof of completeness. Indeed, we 
would have to decide whether the series of particular solutions (17. 3p is the same as 
the obtained general solution. 

Thus, we approach the solution of (17.21) in a different manner. The next theorem 
implies that we can construct a series of particular solutions of (|7.2p using arbitrary 
function z rather than ip, which must be the ratio of two linearly independent 
particular solutions ipi and tp2 of the initial PDE C — Q. 
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Theorem 7.3. Another large class of particular solutions for the PDE 

on I3. Let z E K be arbitrary and non-constant and F ~ F{z) be an arbitrary 
function of z, then 

_ {F{z))x 



q = 



Zy {F{z))y 



and function 



2 

(7.6) /3o(z) = + + 

Zx Z^ IZxZy 

gives particular solutions of (|7.2p . More strongly, if instead of z we use the argu- 
ment F{z), that is l3o{z) is replaced by 

hoiF{z)) , 

we still have particular solutions of i7.2\ l. 

Proof. Can be verified by the substitution. □ 

In other words, on the invariant level we can forget about the fact that t/j must 
be a ratio of two solutions. Of course such trick would not work in general for the 
pre-invariantized system: if we take Darboux Wronskian Formulas and substitute 
arbitrary functions instead of solutions we would not necessarily get a Darboux 
transformation. So the invariantization using gauged evolutions factors out some 
meaningful conditions, which justifies giving them a separate name. 

Using new class of particular solutions discovered in Theorem 17.31 we can find 
the general solution for the invariantized system (see it in Theorem 17. ip in the 
following form: 

Theorem 7.4. Simple Description of all Darboux transformations of bi- 
degree (1, 1). All Darboux transformations of some C = D^Dy + aDx + bDy + c G 
^^[Z)] generated by some M, M = D^ -\-qDy + r G Z^i^*] are parametrized by z Cz K 
an arbitrary non constant and can be written as 



AB 





B 


q = 


"A ' 


h = 


By- A, 


h = 


-Ax + 



where A — and B — — ^ . 

Zx Zy 

Proof. Find z £ K such that q = — — . Then according to (|7.ip we must have 

Zy 

I2 ~ — (ln((7))a;y , which in terms of z has the form given for I2 in the statement of 
the theorem. 

Now we solve (|7.2|) for /a . The solution of this equation can be obtained as the 
sum of particular solution (|7.6p and of the general solution of the corresponding 
homogeneous PDE, 

(7.7) h (— - ^) - h,. + -h.y = . 

\ Zx J 
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Considering I3 in the form I3 = e'^ for suitable J G K[D] this equation can be 
re- written equivalently as 

(7.8) z^_Zy^_.h^J^^^^ 

Zy Zx Zy 

Again we have homogeneous and non-homogeneous parts, and therefore, the 
general solution can be represented as the sum of a particular solution of the non- 
homogeneous part and the general solution of the homogeneous one. 

Using the methods of characteristics one can find the general solution of the 
homogeneous part of (|7.4p . Consider the equality 

dx dy 

-l/Zx l/Zy 

which can be rewritten in the form Zxdx + Zydy = d{z) = 0, and, therefore, z = 
z{x, y) = C, where C is a constant. Therefore, we consider the following change of 
variables: 

e = ^ , 1 
77 = z{x,y) ,J 

which is non-degenerate, since the Jacobian is nonzero: 



Vx Vy 



= rjy = Zy ^ 



d{x,y) 

In the new variables equation (|7.4[) has the form 

and, therefore, the general solution is 

J^H{z) , 

where H{z) is an arbitrary function of z. 

Notice now that since both I^oiz) and Izi){F{z)) are solutions of ()7.2p . therefore, 
their difference is a solution of (|7.7p : subtracting ()7.6|) from /3o(G(z)) we have 

1 (2F'F"' -2,{F"Y\ 
(7-9) hod = ■:^ZxZy \y j = ZxZyG{z) , 

where G{z) is an arbitrary function. Therefore, 

hod = ln(/30d) = Gi{z) + InizxZy) , 

where Gi{z) is an arbitrary function, is a solution of (j7.4p . Therefore, the general 
solution of (|7.4p is 

J ^ Gi{z) + ln{zxZy) + H{z) . 

Correspondingly, the general solution of (17. 7p . which is the homogeneous part of 
the PDE we needed to solve, (|7.2p is / = exp(J), which is 

/ = G{z)zxZy . 

Adding to this the particular solution of the non- homogeneous part, /30, we con- 
clude that 

2F'F"' - 3(F")2 
h = ZxZyG{z) + hoiz) = ZxZy 2{F)^ ^ ^^^^^^ ■ 
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Notice that this expression is exactly l3o{F{z)). Now we proved that 



q = — , 

Zy 



h^h^{F{z)) . 

However, notice that the apphcation of g to z and q to F{z) gives the same result, 
that is q{F{z)) = q{z), the same with I2: l2{F{z)) = hiz). Therefore, we can 
substitute z = F{z) and have the statement of the theorem, which we write in 
terms of z again for convenience. □ 

Comparing Theorem 17.41 with Theorem 17.21 we see that the general solution of 
the invariantized system (see it in Theorem 17. 1|) is much richer than the class of 
the invariantized particular solutions (j7.3p . Does this mean that our hypothesis 
was wrong and that there must be something else besides Darboux transformations 
generated by Darboux Wronskians formulas? In the following theorem lifting our 
results back into pre-invariantized situation we conclude the proof of our hypothesis. 

Theorem 7.5. Given operator C — D^Dy + aDj. + hDy + c G K[D]. Every its 
Darboux transformation generated by A4, A4 = + qDy + r E K[D] is described 
by Darboux Wronskian formulas. 

Proof. Let MC = CiA4 be some Darboux transformations with the gauged evo- 
lution invariants {q, 12,13). According to Theorem 17.41 there exists z & K such 
that 

q - -— - - — 

Zy A 

l2^(^) -(^) =By~A,, 



Cy 



y 



I3 = hoiz) = 



V 
' X 

AB 



2 



where A = , B = — ^ . Comparing this with the invariants given in Theorem l7.2[ 

Zx Zy 

we conclude that there are some Darboux transformation constructed by Darboux 
formulae that has the same gauged evolution invariants. 

Now let us find among the pairs that are constructed using Darboux formulas 
and having these gauged evolution invariants (q, 12,13) those which has the same 
gauge invariants as our initial pair {C, A4). 

Let us denote the gauge invariants of the pair (£, A4) by (mo, ho, qo, Rq). Then 
I2 and I3 can be expressed in terms of (mo, Hq, qo, Rq), see (|6.3p . Therefore, 

2mo = h- {Ro/q)x + Roy , 
2h^ =13 + Rl/{2q) - {Ro/q)x ■ 

Now, since {12,13) are given in terms of z, then invariants {mo, ho) are given in 
terms of z and Rq. 

That is it is enough to find among the pairs that are constructed using Darboux 
Wronskians formulas a pair with the same Rq. 



18 



E. SHEMYAKOVA 



Choose arbitrary functions zi and co and construct an operator C of the form (II. ip 
which has two solutions: zi and zzi and that coefficient c — cq. Using the same 
pair of solutions construct A^' using Darboux formulas, see Theorem 14. 2 1 

Then R' corresponding to the pair (C',M.') can be expressed in terms of z, zi, 
and Co in the form 

R = -— cq + T{z,zi) , 

^2:^1, y ~r ZyZ^ ^ 

where T(z,zi) certain expression depending on z and zi only. This means that for 
every zi we can uniquely find such c that R' = Rq. 

Therefore, for arbitrary zi there exist a Darboux transformation (£', A^') con- 
structed using Darboux formulas for which all gauge invariants (h, k, q, R) of the 
pair are correspondingly the same as those of the initial pair {C,Ai). 

Since those agree, then {£,M) is different from {£',M') by a gauge transfor- 
mation, and therefore, {C,M) can be also constructed using Darboux Wronskian 
formulas. □ 



8. Conclusions 

The present paper closes an essential question for the theory of Darboux transfor- 
mations: Darboux Wronskians formulas are complete for Darboux transformation 
of total order two of operators C = D^Dy + aD^ + bDy + c with non-constant 
coefficients. Since for Darboux transformation of total order there are two famous 
exceptions: Laplace transformations, the case of the total order two has been cru- 
cial. 

We saw that newly introduced transformations of pairs, gauged evolutions may 
have much deeper role than just a tool in proof of our specific problem fTheorem l7.3l 
and the paragraph after it). 

We found a very short invariant description of all possible Darboux transforma- 
tion for C = DxDy+aD.j; + bDy + c generated by in the form Ai = D^ + qDy+r G 
K[D] (Theorem ED). 

Now it is natural to expect completeness of Darboux Wronskians formulas for 
transformations of orders higher than two. We expect this one to be rather difficult 
to prove. Simple repetition and adjustments of the methods and ideas of this work 
would not work. For example, one of the crutial points was the introduction of the 
gauged evolutions, which cannot be defined for pairs {£,M) if Ai has order larger 
than C. 
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